Chaotic scattering through coupled cavities 
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We study the chaotic scattering through an Aharonov-Bohm ring containing two cavities. One of 
the cavities has well-separated resonant levels while the other is chaotic, and is treated by random 
matrix theory. The conductance through the ring is calculated analytically using the supersymmetry 
method and the quantum fluctuation effects are numerically investigated in detail. We find that 
the conductance is determined by the competition between the mean and fluctuation parts. The 
dephasing effect acts on the fluctuation part only. The Breit-Wigner resonant peak is changed to an 
antiresonance by increasing the ratio of the level broadening to the mean level spacing of the random 
cavity, and the asymmetric Fano form turns into a symmetric one. For the orthogonal and symplectic 
ensembles, the period of the Aharonov-Bohm oscillations is half of that for regular systems. The 
conductance distribution function becomes independent of the ensembles at the resonant point, 
which can be understood by the mode-locking mechanism. We also discuss the relation of our 
results to the random walk problem. 

PACS numbers: 05.45.Gg, 73.21. La, 73.23.-b, 05.60.Gg 
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I. INTRODUCTION 

Starting from the study of atomic nuclei, chaotic scat- 
tering has been a topic of intensive research in a large 
variety of systems such as atoms, molecules, quantum 
devices, and microwave cavities [lj. A fundamental ques- 
tion to be asked is how much information is reflected in 
the scattering through random media such as disordered 
and classically chaotic systems. One of the most remark- 
able and promising ideas is to introduce the statistical 
concept into the analysis. The ensemble average over dif- 
ferent realizations of the sample is considered to calculate 
several statistical quantities. A large number of systems 
exhibit universal behavior determined by the symmetry 
of the systems. For this situation, random matrix theory 
(RMT) @,Si has been used to understand the result 
and has played an important role as a standard analytical 
tool. 

Recently, the experimental stage of the chaotic scat- 
tering has been shifted from natural to artificial systems. 
Typical examples are mesoscopic systems [S IS B IS 
such as quantum dots (QDs) and disordered wires. Re- 
cent development of nanotechnology makes it possible to 
fabricate mesoscopic quantum hybrid systems that could 
not be realized before and a lot of interesting interference 
phenomena have been observed under controllable exter- 
nal parameters. Due to the interference of wave func- 
tions, a system made from parts such as the QD, lead, 
and quantum point contact cannot be treated separately. 
Such systems show new interesting phenomena which are 
absent in single isolated systems. Typical experimentally 
fabricated systems are the QD on the Aharonov-Bohm 
(AB) ring [l(|, the side-coupled QD [ll[, and so on. 

The model treated in this paper is two QDs put on the 
two arms of the AB ring. In this so-called "mesoscopic 



double slit system," a lot of interesting phenomena such 
as the AB oscillations and the Fano effects can be ob- 
served by the interference of wave functions transmitting 
through the two arms [HI, [IS 0, HH ■ In the context of 
chaotic scattering, it is interesting to apply the known 
analysis based on RMT [H H3, M, OJ, Sfl, HI HI H| 
to the AB ring system. We study how the interference 
effects appear and the conductance behaves as the func- 
tion of the controllable parameter such as the magnetic 
flux through the ring. 

Our formulation is rather general and the application 
of our result is not limited to the QD systems. It is 
known that microwaves in an irregular shaped cavity be- 
have chaotically and the statistical properties can be de- 
scribed by RMT [13, HI- Based on the formal anal- 
ogy between the Helmholtz and Schrodinger equations, 
the classical waves are simulated as quantum mechanical 
wave functions. Compared with the mesoscopic systems 
in nanoscale, the cavity system is easier to fabricate and 
is ideal for an experimental study. We can also observe 
the Fano effect in this system [26| . 

How can we define the statistical model for the hybrid 
system? For the system of two QDs attached to each 
arm of the AB ring, Gefen et al. [12| considered the case 
when each dot has a single regular level. As a simple but 
nontrivial extension, we treat the case when one of the 
dots has regular levels and the other has random levels. 
RMT is applied to the random dot. 

This model can be viewed as a mixed system of chaotic 
and integrable levels. In single dot systems, such struc- 
ture is employed as an idea to explain anomalous phe- 
nomena such as critical statistics [13] and fractal con- 
ductance [1S|- It is known in the open QD system that 
the several specific levels couple with the lead strongly 
while the other levels couple weakly via strong coupled 



levels [29(. Thus it is too simple to treat the dot as a 
single random matrix and we need to consider the inter- 
nal structure more seriously. Although our model is not 
directly related to such phenomena, it is instructive and 
useful to consider the present ring system as the situation 
where the strong and weak couplings coexist. In this sys- 
tem, the regular transmission in the one arm is affected 
by the random ones in the other arm, and vice versa. 

From a point view of RMT, special attention is paid 
to the universality of the statistical quantities. A natu- 
ral question to be asked in the present model is how the 
universal level correlations described by RMT are mod- 
ified by the regular contribution. Naive expectation is 
that the effect of the regular levels can be safely removed 
by the proper scaling (unfolding) [3(j. It is known that 
the effective theory is written in terms not of the mi- 
croscopic parameters but of the transmission coefficients 
[ill. However, in the present system, the effect is am- 
plified by multiple scatterings through the ring and gives 
highly nontrivial results. 

Now that our model has been described, we must refer 
to the work by Clerk et al. [311]. They considered many 
resonant levels in a single dot and RMT was employed 
for their distribution. The regular component to the S 
matrix expressing the direct nonresonant path through 
the dot was used to find the Fano resonances. For each 
resonance, the Fano parameter q was calculated and the 
statistical distribution of q was defined over the reso- 
nances. On the other hand, in our case, only the single 
resonant level is present regularly and it is affected by 
random levels. Thus our attention is fixed on the single 
regular resonance. To discuss the statistical properties 
of the transport we must prepare different realizations 
of the random dot. The ensemble average is defined in 
terms of such realizations. 

The outline of this paper is the following. The AB ring 
model is defined in Sec|TTJ We define the random Hamil- 
tonian model in Sec lII Al A model based on the random 
S matrix is also defined in Sec lII Bl and the relation to the 
random Hamiltonian model is discussed. In Sec lIIII we 
calculate the average of the S matrix based on the ran- 
dom Hamiltonian model. As a result the mean part of 
the conductance is calculated. It is not enough to calcu- 
late the conductance including the quantum fluctuation 
effect and we develop the supersymmetry method [32J in 
Sec II V Al to calculate the full conductance. The results 
of the conductance are shown in Sec lIVBl We also study 
the AB oscillations in Sec lIV Cl and the Fano effect in 
Sec llVDl The fluctuation effects can be best seen in the 
conductance distribution functions, which are studied in 
Sec|V] Since realistic situations are not ideal and phase 
breaking effect is present [25, 33], it is important to con- 
sider the dephasing effect theoretically. We consider it in 
Sec lVIl using a simple imaginary-potential model. Sec- 
tion [VlT] is devoted to discussions and conclusions. Part 
of the results were published in a preliminary report [3~ij . 



dot1 



Lead 




Lead R 



dot 2 

FIG. 1: Schematic drawing of our model. Dot 1 with a reso- 
nant level Ei and dot 2 with a random Hamiltonian H2 are 
connected by leads, w denotes a dot-lead coupling matrix. 



II. MODEL 

A. Random Hamiltonian approach 

We consider the AB ring system depicted in FigJTJ The 
upper dot (dot 1) has a single resonant level, and the 
lower dot (dot 2) has random levels and is treated by 
RMT. It is known from scattering th eory that the S ma- 
trix of the system is written as [g, 0, [H, HH 
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where H denotes the Hamiltonian matrix for dots and w 
the dot-lead coupling matrix. H can be written as 



H = 



Ex 
H 2 



(2) 



where E\ is the fixed energy level for the dot 1 and H2 
is the random Hamiltonian for the dot 2. The size of 
H2, N, is taken to be infinity to find the universal re- 
sult. We note that the total size of H is 1 + N, It is a 
straightforward task to extend the size of the upper dot 
Hamiltonian to arbitrary values and here we consider the 
minimal size 1. As another simplification, we consider 
the 2x2 (unitary) matrix S, which means that the left 
and right leads have a single channel, respectively. It 
is believed that the quantum interference effect becomes 
maximal in this case [35| . Then the dot-lead coupling 
matrix w is the (1 + N) x 2 matrix and can be written as 
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where (1L) refers to the coupling between the dot 1 and 
lead L, and so on. 

The conductance measures the transmission from the 
left to right lead and is defined by 0, 0, EE] 



5=<|Si 2 | 2 >, 



(4) 



where ( ) denotes the ensemble averaging of the random 
Hamiltonian H 2 . We employ the Gaussian ensemble [3j 
and the probability density is given by 



P(H) = Cexp 



2NA- 



(5) 



where A is the mean level spacing, and C is the normal- 
ization constant. In the following calculations, we mainly 
consider unitary symmetry, which means that H 2 is Her- 
mitian and no additional condition is imposed. 

The result of the conductance depends on the choice 
of the dot-lead coupling w. Although this matrix w has 
4N degrees of freedom, there is no need to specify them 
completely. After the averaging, the effective degrees of 
freedom becomes finite. Generally, it is 6 and we restrict 
our discussion to the special case of 4 (see below). 



B. Random S matrix approach 

Equation (fTJ) is a useful formula to relate the Hamil- 
tonian to the S matrix and can be used for the present 
coupled system. It is convenient to express the S matrix 
in terms of the K matrix defined by 



S = 



1-iK 
1 + iK' 



(6) 



K is expressed as the sum of contributions from dot 1 
and 2: 



1 



K = K 1 +K 2 

K = 7TuA 

K x = 
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This simple relation implies the sum rule of the S matrix 
1-5 1 - Si 1-S 2 



1 + S 1 + Si 



1 



(8) 



where Si (S2) is the S matrix for dot 1 (2). It is instruc- 
tive and useful in the following numerical calculations to 



derive the explicit representation using the matrix ele- 
ments. Defining each S matrix elements as 



S = 



r t' 
t r' 



we obtain, for example, 



h(l + s 2 + r 2 + r' 2 ) + t 2 (l + si + n + r[) 

9 + 3(n + r[ +r 2 + r' 2 ) + s : + s 2 
3(rir 2 + r[r' 2 ) + {r x r 2 + r 2 r[) - 4(*i4 + t 2 t[) 



-(ri + r[)s 2 - (r 2 + r' 2 )s 1 + sis 2 



(10) 



where Sj = det Si — — (i = 1,2). Thus the total 
transmission t is not equal to ti+t 2 , rather including non- 
linear effects due to multiple scattering through the ring. 
Such multiple scattering effects are put together with in- 
terference due to random scattering and give highly non- 
trivial results for the conductance g = {\t\ 2 ). 

Another way of representing the total S matrix is to 
separate the S matrix of the system into the upper and 
lower dot parts and the left and right fork parts [12l . [l3| . 
Choosing the fork matrices in a proper way, we can find 
the same expression of t as in Eq. (fT0| . 

The conductance can be calculated by taking the en- 
semble average over S2 determined by the random Hamil- 
tonian H 2 . Instead of doing that, we may disregard the 
detailed structure of 5*2 and impose randomness directly 
on S 2 , simulated by the circular ensembles It is well 
known that the random S matrix approach is equivalent 
to the random Hamiltonian approach if we use the Pois- 
son kernel LL7l 



P (S)d^(S) cx 



|det(l-5(5)t)| 
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dnf,(S), (11) 



where dfip(S) denotes the invariant measure for the S ma- 
trix and is used as the measure for the circular ensemble. 
j3 is the index for the universality class. (3 = 1,2, and 4 
for the unitary, orthogonal, and symplectic case, respec- 
tively. (S) is the averaged value of S which is treated 
as an input parameter and is determined by the random 
Hamiltonian model. The total S matrix is constructed 
by the sum rule ((5|) and the conductance is expressed by 
|t| 2 where t is given by Eq. tfTOf . By taking the circu- 
lar ensemble average with the weight Pp(S 2 ), we obtain 
the conductance g which is the same as that obtained 
by the random Hamiltonian approach. The equivalence 
of both approaches was shown in Ref . [22| . The random 
S matrix approach has a great advantage for numerical 
calculations because there is no need to take the ther- 
modynamic limit N — > 00 and one may consider 2x2 
random matrices S^. 

Alternatively, we can parametrize the S matrix in 
terms of the K matrix ((7]). Then the expression of the 
conductance becomes much simpler than Eq. (|10p as we 
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show in Sec|Vl The disadvantage of this parametrization 
is that the K matrix is Hermitian and the matrix elements 
are not compact, which is inconvenient for the numerical 
calculation. Thus we employ the S matrix parametriza- 
tion (|10p with compact variables for most of the numer- 
ical calculations. 



III. AVERAGED S MATRIX 

As we have shown in Eq. ([7]) , the K matrix is written as 
the sum of the regular (dot 1) and random (dot 2) parts. 
Thus, to get the averaged K matrix, we may consider the 
ensemble averaging of the random part. We know from 
RMT that the averaged Green function for the Gaussian 
unitary ensemble is given by [3j 



where 



H, 



cosz 



NA 



TlE 

2NA' 



(12) 



(13) 



N is taken to be infinity while E/Ais kept finite. Then 
we have e~ tz — > — i and the averaged K matrix is given 
by 



(K) = 



1 



Ei 



7i 



NA 



72, 



(14) 



where j{ = ttw^w^ (i — 1,2). It is important to note 
that the result depends on the dot-lead couplings u/ 1,2 ' 
through 7i ; 2. 

For the regular dot, the most general form of 71 is 



71 



TTW 



(1Z) 



nw (lR)* w (lL) 



ly(lfl) 

{1R)* W (1R) 



\/rii?Ti L e 



VrizTifle' 

r 1R 



(15) 



where Til (r^) turns out to be the level width for the 
left (right) coupling of the dot to the lead and is a 
phase. We assume the symmetric coupling Til = Vm 
for simplicity and use 



where 



71 = 



2 V e 1 



(16) 



(17) 



This matrix satisfies $ 2 = $ and is diagonalizcd as $ — > 
diag(0,l). 

On the other hand, for the random dot, the form of 72 
is slightly complicated. It is written as 



72 



,,(2L)t 7 .,(2L) 



,,(2i)t,.,(2fi) 



7rw (2fl)t w (2L) nw (2RK w (2R) 



(18) 



and ti/ 2fl ) 



Since w {2L) 
see that the 

| w (2L)t w (2i?.)|| w (2fl) t , 



are Af x 1 matrices, we 
relation | w (2i)t w (2L) | | w (2fl)t w (2fl) | > 

\ W \ 2L )\ holds. The equal sign holds 
when vj( 2L ) — w^ 2R ^ or N = 1, the latter is the case 
for u/ 1 ) . Thus we need the additional parameter for the 
parametrization of 72. Assuming the symmetry of the 
left and right coupling again, we obtain the form with 
the level width T2 as 



AT 2 



72 



1 ae 

D -i4> x 



i<f> 



(19) 



The parameter a reflects the above mentioned inequality 
and < a < 1. We note that the same phase <j> appears 
in 71 and 72, but the sign is opposite to each other. This 
phase affects the transmission part of the S matrix and 
can be identified with the AB flux through the ring pjl 

Da. 

Using this parametrization, we can write 



(K) 



1 



iX 



1 
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1 
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where 



E-E x v _TrT 2 
Ti ' A 



(20) 



(21) 



The energy e represents the distance from the resonance 
point and X is the ratio of the level width to the mean 
level spacing of the dot 2. Thus this model is described 
by four parameters e, A, a, and <p. 

For the random dot, the elements of the dot-lead 
coupling distribute randomly and the summation 



E 



N (2L)*(2R) 



can be small when the random phases 



of w( 2L ^ and u/ 2i? ) almost cancel out. This means a is 
vanishingly small. On the other hand, the summation 
can be finite when the left and right dot-lead couplings 
are correlated mutually. This results in direct nonreso- 
nant reaction We first discuss the case of a = for 
simplicity. The averaged K matrix takes the form 



1 iX 
_$ _ 

e 2 



(22) 



The finite-a effect is discussed afterwards. 

Now we go back to the S matrix. The averaged S ma- 
trix is simply obtained by using the averaged K matrix, 



(S) = 



l + (K) 



1 - 



x 



2i 



(23) 



1 + | (l + f)[(l + f)e + <] 

This is justified by the saddle-point analysis of the non- 
linear sigma model described below. We define go = 
\(Si 2 )\ 2 , which is the conductance if we can disregard 
the quantum fluctuations. It is given by 



1 



9o 



2 , 



(£-£i) 2 



1+4 



(24) 



5 



The result shows that the level width Ti for the dot 1 and 
the conductance are renormalized by the factor 1/(1 + 
X/2). 

For later use, we define the transmission coefficients as 
T = 1— (5)(5) f 

2X 2X 



1 



x \ 



(1 + f) {E-E 1 Y + (^ K 



(25) 



This matrix can be diagonalized to find the eigenvalues 



Ti 



2X 



1 + 4 



To 



2X 



(1 + f) 



l 

7? 



(26) 



Note that < T 2 < T x , T 2 = T x at \E - E x \ -> oo, and 
T 2 = at E = Ei. At X = 2, T x takes the maximum 
value, T x = 1, and the transmission through the random 
dot becomes ideal. 

In conclusion of this section, we found the averaged S 
matrix ([23]) and the conductance ([24]) . Of course, this is 
not the final result of the averaged conductance. We just 
calculated the mean part go = |(<Si2)| 2 which is different 
from the original definition We must examine the 
fluctuation part 5g = g — g — (\Si2\ 2 ) — |(<Si2)| 2 - 



IV. CONDUCTANCE 



This is easily solved with the proper boundary condition 



where we took the limit N — > oo keeping E/A finite. As 
a general solution including the saddle-point manifold, 
we can write 



Q = -ia, a = VAV, 



(31) 



where V is the 4x4 supermatrix and satisfies VV = 1. 
The symmetry of V is determined in the standard way 

Now the generating function reads 



(Z) = 
F = 



Vae 



str 



4(1+JV) 



In 



str h In 



1 + iA 



1 171 



E-Ei 



7i 



iVA 



72 (T 



(32) 



where we assumed 71 = Ti$ and 72 = NT2/2. Since 
the matrix sizes of a and 7 are 4 and 2, respectively, the 
total size of the superspace in the last expression is 8. 
We finally obtain 



A. Supersymmetry method 

We derive the nonlinear sigma model for the coupled 
system to calculate the fluctuation part of the conduc- 
tance. According to the supersymmetry method (T5 , |32| , 
the generating function for the product of Green func- 



tions G(E) 
by 



l/(E-H 



iirww 



) and G\E) is defined 



T>(ip,ip)exp iip(^E + iAirww^ - H^jip , (27) 



where ip has 4(1 + N) components coming from su- 
persymmetry (bosons and fermions), retarded- advanced 
structure, and Hamiltonian space. A = diag(l,— 1) in 
retarded-advanced space. Following the standard proce- 
dure, we introduce the Hubbard-Stratonovitch field Q to 
write the averaged generating function as 



(Z) 



VQ exp< —str 



-4(1+N) 



hi 



E + iAirwiv' 









7T ^ 



- 2 Str 4<9 2 



(28) 



where Q is a 4x4 supermatrix. "str" denotes supertrace 
and the subscript indicates the size of superspace. When 
w = 0, the saddle-point equation is written down as 



iVA 



1 
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NA 



Q 
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The nonlinear sigma model (|33[) can be written in terms 
of the transmission matrix T (|25p and the microscopic 
fundamental parameter X does not appear in the expres- 
sion explicitly. This is a manifestation of the universality 

Equation (|33[) is for systems with unitary symmetry. 
In the same way we can derive the nonlinear sigma mod- 
els for the orthogonal and symplectic symmetry classes. 
Then a becomes an 8 x 8 supermatrix and the additional 
symmetry due to time-reversal invariance is imposed [32| . 



B. Conductance 

In the nonlinear sigma model approach, the averaged 
conductance is calculated by performing the integration 
of the a matrix. The mean part of the conductance go in 
Eq. is easily obtained by neglecting the fluctuation of 
the a matrix as a = A. To find the fluctuation part of the 
conductance Sg = (IS12I 2 ) — |(<Si2)| 2 , we must take into 
account the contribution from the saddle-point manifold 
parametrized by the V matrix. This calculation is highly 



FIG. 2: Conductance vsl = 71T2/A. The thick (thin) lines 
are analytical results of the total conductance g (mean part 
go)- Inset : Comparison of the analytical (denoted by lines) 
and numerical (dots) results for the fluctuation part Sg = 

9 ~ go- 



complicated and we refer to the Appendix „ for details. 
We finally obtain 



Sg 




where are the eigenvalues of the transmission matrix 
T given by Eq.(|26j). 

We first examine the two limiting cases, |e| — > 00 and 
e = 0. The limit |e| — > 00 means that dot 1 is detached 
from the system and the S matrix is given by Si = 1. In 
this case, T x = T 2 = 2X/(1 + X/2) 2 and we recover the 
known result 36] 



(35) 



In the other limit e = (E = Ex) the energy coin- 
cides with the level in dot 1 and the perfect transmission 
through dot 1 is achieved. Then T 2 = and we obtain 



69 = T . 



(36) 



We see that Eq. (|3"5"|) is larger than Eq. (f3"6"|) , which means 
that the fluctuation effects are reduced as we approach 
the resonant point. For intermediate values of e, Eq. (f3~4")) 
cannot be written in terms of T1.2 only in contrast to 
Eqs. (j3"5)l and (f3"6")> . This is because the source term to 
calculate the conductance depends on e and X explicitly, 
although the nonlinear sigma model itself can be written 
in terms of T, as shown in the Appendix „ 

These results are checked by numerical calculations. 
We use the formula (|10p for the transmission matrix. Si 
is given by Si = (1 — iK\)/(l + iK{) with K x = <f>/e, and 



FIG. 3: Conductance vs e = (E - Ex)/Yi- The thick (thin) 
lines are the total conductance g (the mean part go). Inset : 
Fluctuation part Sg = g — go- 
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FIG. 4: Numerical results of the conductance g(e) for orthog- 
onal and symplectic ensembles. The result for the symplectic 
case is normalized to unity. 



the random S matrix S2 is treated statistically by using 
the Poisson kernel (fTT|) . We take the ensemble average 
over more than 10 6 samples of the S matrix. 

In Fig. El X dependence of the conductance is shown 
for several values of e. go shows a peak at X = while 
Sg takes a maximum at X = 2 as shown by the thin lines 
and the inset in Fig. [2 respectively. As e — > 00 go {Sg) 
is monotonically decreasing (increasing) and the result 
rapidly approaches Eq. (|35[) . The numerical result agrees 
with Eq. (|34p in a highly accurate way, which shows the 
equivalence of the random Hamiltonian and random S 
matrix approach. 

e dependence of the conductance is shown in Fig. [3] 
A resonance peak appears at e = 0, reflecting transport 
through the regular dot 1. This peak structure, how- 
ever, changes qualitatively as a function of X. For small 
X the peak is convex and the peak height decreases on 
increasing X. When X = 2, g is independent of e. In- 
creasing X further, we find that the peak turns into an 
antiresonance and g decreases monotonically. The result 
for X — 2 corresponds to that of the circular unitary en- 
semble because (S2) = 0, and the Poisson kernel Pp{S2) 
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FIG. 5: Conductance vs <f> for X = 2. The lower figure is 
for the orthogonal ensemble and upper for symplectic. No 
oscillations are observed for the unitary case. 



becomes unity. As we see in the inset of FigEl $9 a t 
the resonant point is relatively small and the quantum 
fluctuation effect smooths the resonance. 

For comparison we calculate g as a function of e for the 
orthogonal and symplectic ensembles numerically For 
the orthogonal case, the Hamiltonian has time-reversal 
invariance and the matrix elements are real. For symplec- 
tic, the Hamiltonian becomes a quaternion real matrix 
The results are shown in FigQ] When < X < 2, the 
resonance is enhanced (reduced) for the orthogonal (sym- 
plectic) ensembles. At X = 2, the orthogonal ensemble 
gives a resonance while the symplectic ensemble gives an 
antiresonance. When X — 10, we see that antiresonancc 
is reduced (enhanced) for the orthogonal (symplectic) en- 
semble in contrast to the case of A < 2. 

Away from the resonance, the quantum fluctuation ef- 
fect becomes larger as the number of degrees of freedom 
of random variables increases. We note that the number 
of degrees becomes maximum when (3 = 4 and minimum 
when f3 = 1. We can also see that the conductance at 
the resonant point is independent of the choice of the 
ensemble. This result is discussed in detail in SecfVl 



Aharonov— Bohm oscillations 



In systems with unitary symmetry, since the scatter- 
ing in the random dot randomizes the phase of the am- 
plitude, the result becomes independent of the AB phase 
(\>. This is not the case for the orthogonal and symplectic 
systems and the oscillations can be observed. However, 
the period of the oscillation is different from that for regu- 
lar systems. This can be understood from the expression 
of the transmission t in Eq. fTUl) . The phase is included 
in that expression as 



At x e 



i(j> 



Bt 2 e 



-i<f> 



(37) 



where A, B, C, and D are phase independent contri- 
butions. If we neglect the multiple scattering effect the 
total transmission is approximated by t ~ t\e % ^ + t 2 e~ 1 ^. 
Then the conductance is given by 



g ~ \he^+t 2 e-^\ 2 
= \h\ 2 + \t 2 \ 2 + t l t* 2 e 214 



t\t 2 e 



-2i<j> 



(38) 



We see that the third and fourth terms of the right hand 
side give oscillations with the period tt. However, these 
terms vanish after the random averaging. The contribu- 
tions going around the ring twice give oscillations with 
the period 7r/2 and survive after the averaging. Such con- 
tributions come from expanding the denominator. Thus, 
in the orthogonal and symplectic systems, g depends on 
the AB phase due to the multiple scattering inside the 
ring. The period of the AB oscillations becomes half of 
that for the regular systems. This effect can be inter- 
preted as a kind of the Altshuler-Aronov-Spivak effect 
^37'] for cylinder systems. In a ring system, it was dis- 
cussed in Ref.[38| that the period of the oscillation be- 
comes half a flux quantum by the self-averaging effect. 

We show the numerical results in Figj5] for the orthog- 
onal and symplectic ensembles. The period of the oscil- 
lations is 7r/2 as we discussed above and the difference 
between these two results is that the conductance be- 
comes minimum (maximum) for orthogonal (symplectic) 
at (j) = 0. This can be understood by the standard mech- 
anism of weak localization 13911 . 



For regular ring systems, it is well known that the AB 
oscillations are observed by applying the magnetic flux 
through the ring. Since the flux is a tunable parameter 
it is an important method to control the system. Our 
interest is how the effect of the AB flux can be observed 
in the present random system. Can the AB oscillations 
survive after the random averaging? 



D. Fano effect 

The Fano effect is induced by the correlation of the 
resonant and direct path [3l], [4(| ■ The direct path can 
be described by the parameter a in Eq. (fI9|) . If we keep 
this parameter in Eq. (|20p , the averaged S matrix is given 
by ' 



(S) 



(l+l±aX)(l + i^X) + (l + 



1 — a cos 2c 



x 



8 



5^ - ale-'* 1 - i^-X 2 - 



r 3 1 — a cos 20 yi 

(39) 
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in <5g as shown in the inset of Fig [31 the asymmetry is 
opposite to that of go- As a result, the total conductance 
becomes symmetric. This result means that the Fano 
effect appears not on g but on the mean part go and 
the fluctuation part Sg, respectively. We note that the 
asymmetric form is obtained when Re q ^ 0. The effect of 
the imaginary part of q keeps the conductance symmetric. 
We can conclude that the real part of the Fano parameter 
does not affect the total conductance. We confirmed that 
the symmetric conductance is obtained for the orthogonal 
and symplectic classes as well. 



FIG. 6: Conductance vs e = (E - E\)/T\ at a = 0.7 and 
4> = — 7r/8. The thick (thin) lines are the total conductance g 
(the mean part go). Inset : Fluctuation part Sg — g — go- The 
total conductance g is obtained numerically and the mean 
part go is plotted by using Eq. (|40)) 



The mean part of the conductance is derived from this 
expression as 



a 2 X 2 



.90 = 



(1 + l±a X f (l + ^X) 2 



|e - ei + qTi 



(E - Ei) 2 



-, (40) 



(x + i±±x)*(i+^x) 
where q is the Fano parameter 



-F 2 
1 l 



aX 



(41) 



Thus the Fano effect appears when a/0. The additional 
condition <f> ^ is required to obtain a finite real part of 
q. Then the asymmetric conductance form is obtained. 
At the limit \E—E\\ — > oo, go has a finite contribution in 
contrast with Eq. (j2"4|) . This means that there is a direct 
regular coupling between the left and right leads through 
the random dot. 

This Fano effect also appears on 5g. The numerical 
result in Figj6] shows that the Fano parameter for 5g is 
the same as Eq. (l4"TT) . Since the antiresonance appears 



V. CONDUCTANCE DISTRIBUTION 
FUNCTIONS — MODE-LOCKING EFFECT 

In the previous section we focused on the averaged con- 
ductance. It is well known that disordered systems show 
strong fluctuation effects, which mean that the square of 
the conductance and even the higher moments become 
relevant to characterize the system. To discuss the effects 
of the fluctuations, here we calculate the conductance 
distribution function P(g) — (S(g — |Si2| 2 ))- We show 
the analytical results when e = which show universal- 
ity among the ensembles. We also report the numerical 
results. 

The expression of the conductance distribution be- 
comes simpler if we use the K matrix representation as 
we mentioned in Sec lII Bl K is a Hermite matrix and 
K = Ki + K 2 with 



1 



i 



—i(j> 

i 



Kn 



ai 
tie-' 



be^ 

2 



(42) 



ai.2 are real, and b depends on the universality class and 
is expressed as 



b = 



b 



bo 

b + ib\ 
b\e\ + 6 2 e 2 



for (3 
for f3 
b 3 e 3 for 13 



(43) 



where 60,1,2,3 are rea l an d ei.2,3 quaternion matrices de- 
fined by ej — iuj with the Pauli matrix <jj (j = 1,2,3). 
The conductance is expressed in this paramctrization as 



l+4& O €cos20+4|fc| 2 e 2 



[60 cos 20- 



l+(l-aio 2 + |6| 2 )e +[l+(a 1 +a 2 )£] 2 



l+4(bpe cos 20-bi sin 20)+4|&|' 



% cos 2(f) — bi sin 20- 



i^a + (l-aia2 + |f>| 2 )«r + [l + (ai+a2)e] 2 



tr 



[(fro cos 20— r<T3 sin 2</ 



+ (l-a 1 a 2 + \b\ 2 )e]' + [l+(a 1 +a 2 )ef 



for (3=1, 
for (3 = 2, 
for (3 = 4, 



(44) 
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where \b\ 2 = an d 7,2 — b\ + b\ + b\. We note that the conductance for (3 — 4 is normalized to unity. This 

expression is averaged by the generalized circular ensemble (Poisson kernel) 

t ! >| (/3+2)/2 0-X 

Pfi(3 2 )dnp(Su) oc I — — ——2 I d ai da 2 Y[ db i: (45) 

I [f + x(- a i«2 + H 2 )] +(«i + a 2 ) 2 J f=o 



where we used (S 2 ) = (1 - X/2)/(l + X/2). The numeri- 
cal results using the Metropolis algorithm j4l| are shown 
in FigEat = 0. 

The results at large e are interpreted as the single ran- 
dom dot case. This case was discussed in Ref.[18j using 
the random Hamiltonian approach and the analytical re- 
sult for the unitary system was obtained. In the random 
S matrix approach, the case of the perfect transmission 
X = 2 was obtained in Ref.|20] as 



P(9) = f 9- 1+0/ \ 



(46) 



and other cases were discussed in Ref.[21|. The case of 
e = 10 is enough to find a large-e result and we find a 
good agreement with their results. 

In the case of e = 1, we can clearly see how the random 
dot significantly affects the distribution function. If we 
increase X, a single peak at small X turns into a broad 
one and a different peak around g — is formed at large 
X. 

It is interesting to see the results at e = which are in- 
dependent of the choice of the ensemble. In this case, the 
conductance distribution can be calculated analytically. 
The conductance is written as 



1 



1+ (bo 



(47) 



and the distribution function is obtained from the ex- 
pression 

0-1 



P(g) = C I d ai da 2 Y[ dbi 
J 1=0 

xS[g 



(/3+2)/2 



[f + x(- a i a 2 + H 2 )] +{ai+a 2 ) 2 



(48) 



where C is the normalization constant. We perform the 
integrals and obtain 



P{g) 



1 



1 



7r^/g(l-g)i(l-g) + ±g 



(49) 



This result agrees with the numerical ones in FigtZJ The 
reason why this result becomes independent of (3 can be 



considered as follows. In Eq. (|42[) , all the matrix elements 
of K\ are divergent at a = 0. When <\> = 0, this diverging 
term belongs to the member of the orthogonal ensemble 
and affects the variables ai j2 and 60 in the second term 
K 2 which are common to all the ensembles. Then the 
effective modes are locked on those for the orthogonal 
class and the conductance (|T7|) becomes independent of 
the rest of the parameters b\ 2 3. 



Equation (|49p with X = 2 appears in the problem of 
the classical random walk [i^] and is known as the arc- 
sine law. Consider the one-dimensional classical random 
walk starting at the origin. The walker can move to ei- 
ther one of its two nearest neighbor sites with the equal 
probability p = 0.5. After the TV-step walk, we count the 
number of the events which the walker was in the posi- 
tive axis M. Then the distribution function of g = M/N 
approaches Eq. (|4"9")) with X = 2 as N — > 00. It can 
be considered that the walker at the positive (negative) 
direction corresponds to the transmission (reflection) to 
the left (right) lead in our model. Due to the presence of 
the resonant path through the dot 1 , a particle transmit- 
ted through the dot 2 can go to either left or right lead 
with equal probability. The particle reflected by the dot 
2 can go either way as well. Thus the particle entered 
from a lead forgets where it came from. Such a process 
can be interpreted as a random-walk-like one and gives 
the same distribution function. It is interesting that the 
asymmetric random walk with the probability p ^ 0.5 
can be described by our model with 1/2. Since the 
analytic form is not known in the asymmetric random 
walk, our result may be useful for understanding the re- 
sult. It is also known that the same distribution function 
appears in the problem of the continuous-time quantum 
walk 

When the phase </> is finite, K\ does not belong to the 
member of the orthogonal ensemble and the results can 
depend on the choice of the ensembles. We numerically 
found that the orthogonal and unitary cases are indepen- 
dent of (j> and the result (|49[) is kept unchanged. For the 
symplectic case, we found that the result depends on <fi 
and Eq. (|49p does not maintain anymore. The numerical 
result for X = 2 and e = is shown in Fig[8l Remarkably, 
all plotted curves give the averaged conductance g = 0.5 
and the phase dependence appears only for the conduc- 
tance fluctuations. We also see that plotted curves has 
a nonanalytic point at around g = 0.5, which implies a 
nontrivial mechanism due to the phase coherent effect. It 
is not clear how this happens and further study is needed 
to clarify the underlying mechanism. 
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FIG. 7: Ensemble dependence of the conductance distribution 
functions at cj> = 0. The curves at e — are well fitted by the 
analytical result (|49[) 



VI. DEPHASING 

In the Hamiltonian approach, the dephasing effect can 
be modeled by introducing the imaginary part to the en- 
ergy 
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2t' 



(50) 
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FIG. 8: Conductance distribution function of the symplectic 
system for several values of <j>. We set X — 2 and e = 0. 
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FIG. 9: Conductance distribution function for unitary system 
for several values of p = 1/Ar. 



This method is equivalent with that of Ref. [l8[ where the 
imaginary part of the Hamiltonian was introduced. In 
the supersymmetry method, this effect can be described 
by the additional term of the sigma model [l|[ 



F T = — strcrA. 
At 



(51) 



This term makes the massless "diffusion" modes on the 
saddle-point manifold massive and reduces the quantum 
fluctuations. See the Appendix „ for details. 

It is well known in the S matrix approach that the de- 
phasing effect can be described by the Biittiker's voltage 
probe model [44|. A fictitious voltage probe eliminat- 
ing the phase coherence is attached to the dot and is 
described by an enlarged S matrix. 

Brouwer and Beenakker showed that the voltage-probe 
model at a certain limit becomes equivalent to the 
imaginary-potential model and found the modified Pois- 
son kernel in the random S matrix approach [351 ] . Here 
we investigate this limit using the imaginary-potential 
model. Since the dephasing effect to the regular dot is 
trivial, we include the effect in the random dot only. 

In FigOU the numerical results of the conductance dis- 
tribution function using the random Hamiltonian model 
are shown. We add the dephasing term, p = 1/Ar with 
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the phenomenological dephasing rate r, to the Hamilto- 
nian. The matrix elements of the dot-lead coupling 
are chosen randomly so that there is no direct nonreso- 
nant reaction a = 0. The size of the random Hamiltonian 
is taken 10 2 and the averaging over 10 5 samples is carried 
out. 

As p increases, the curve transforms into a single peak 
structure. The peak point corresponds to the mean part 
of the conductance go, which is close to zero for the upper 
graph and 0.25 for the lower one. We can conclude that 
the dephasing effect only affects the fluctuation part. We 
also confirmed that our numerical result based on the 
random Hamiltonian agrees with that of the random S 
matrix model in Rcf. 35]. 



VII. CONCLUSIONS 

We have discussed an AB ring system with regular and 
random cavities. We found that the quantum fluctua- 
tion effect plays an important and crucial role and sig- 
nificantly affects the conductance. The main results are 
summarized as follows: (i) The averaged conductance is 
divided into two parts. The mean part has the Brcit- 
Wigner resonant form renormalized by random effects. 
The quantum fluctuation part has an antiresonance form 
where the quantum effects become minimal at the reso- 
nant point, (ii) For the orthogonal and symplectic en- 
sembles, the AB oscillations are found and the period 
of the oscillations is half a flux quantum. The positive 
(negative) magnetoconductance are obtained for the or- 
thogonal (symplectic) ensemble because of the multiple 
reflections inside of the ring, (iii) Depending on the pa- 
rameter choice, the Fano effect can be observed. This 
effect appears in the mean and fluctuation parts, respec- 
tively, and a symmetric form is obtained for the total con- 
ductance, (iv) The conductance distribution functions 
clearly show the influence of strong fluctuations. The 
distribution function at the resonant point, Eq. (|49p . does 
not depend on the choice of the ensemble, which can be 
understood by the mode-locking mechanism. The form 
of the distribution function implies a relation to the ran- 
dom walk problem, (v) The dephasing effect simulated 
by the imaginary-potential model reduces the fluctuation 
part only. 

The result of the averaged conductance in Fig[3] shows 
that the total conductance as a function of the energy 
takes a broad distribution. The form of the total con- 
ductance is determined by the competition between the 
mean ([24]) and fluctuation (|34| parts. At large X we can 
observe the antiresonance. 

Separating the mean and fluctuation parts is crucial to 
understand the obtained result. For example the Fano ef- 
fect is found in both parts, while the total conductance, 
the sum of them, becomes symmetric. We also found 
that the dephasing effect suppresses the fluctuation part, 
which means that the cancellation is incomplete and the 
asymmetric form can be obtained in a system with de- 



phasing. 

The most striking result can be seen in the calculation 
of the conductance distribution function. At the resonant 
point, the effective modes of the K matrix in Eq. ([42|) are 
locked to those for the orthogonal ensemble. Only the or- 
thogonal modes are amplified by the multiple scattering 
through the ring and we can find the ensemble-insensitive 
result. This result suggests a possibility of controlling the 
ensemble dependence of random systems by the resonant 
singularity embedded in the systems. 

Our results for a coupled system show that the non- 
trivial phenomena which are absent in the single system 
can be observed in the hybrid system, which opens a 
new direction for theoretical and experimental studies 
of chaotic scattering. In this paper we only considered 
the regular-random coupled system. It is interesting to 
see more complicated systems such as a random-random 
system and triple coupled cavities, and so on. A study 
of the series coupled random dot can be seen, e.g., in 
Ref . pH . To the best of our knowledge, there is no sys- 
tematic study on the parallel coupled system. It will be 
discussed in detail in a future publication [461 ]. 
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APPENDIX: CALCULATION OF THE 
CONDUCTANCE 



We calculate the conductance using the nonlinear 
sigma model with unitary symmetry, Eq. (|33p . The first 
step to do is to represent the conductance as an inte- 
gral of the er matrix. This is the standard prescription 
discussed in detail in Ref. [l6[ and we have 



|S 12 | 2 > = (str lk±±±K 



1 



1 + iAK 



, , 1 - A r> 1 
xstr k K- 



l + iAKj 21 / F 



1 + iAK J u 



X | fc 1 -— if- 



1 + iAKJ 22/ F 



(A.l) 



where k = diag(l,— 1) in superspace, ( }p denotes the 
integration over a with the weight e~ F , and 



E — E\ NA ' 

1 iX 

= -<E> (7. 

e 2 



In the second line we used a — 0. 



(A.2) 
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The second step is to parametrize the supermatrix a. 
We use [HI 

a = U anU. 



= 
U = 



COS V 

-i sin/ 



% sm v 
— cos ( 



R A 
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RA 



where 



z0 B 

o F 



(A.3) 



(A.4) 



Bl- 



and the integration range is given by < 9 b < oo and 
< 9p < ir. U includes the anticommuting Grassmann 
variables and can be written as 



U = U\U2, 



u\ = exp 



irj 
-ir]* 



BF 



»2 



v = exp 



e tpi o 
e^ 2 

K 



-K* 



BF 



BF 



(A.5) 



where 77 and \ are Grassmann variables and the range of 
the real variables (pi t 2 is given by < ipt t 2 < 2n. The 
invariant measure of this parametrization is 



Da = Cd9Bd9Fd(pidip2dr)drf dudn* 

sinh 9 b sin 9 F 
x 

(coshes — cos# F ) 2 ' 



(A.6) 



where C is the normalization constant. In this 
parametrization, we can write 

e _ F = [1 - - cosfl F )] [1 - _ cosg F )] 

[l + (cosh 9b - 1)] [1 + ^(cosh6» B - 1)] ' 

(A.7) 

The last step is to carry out the integrations. This cal- 
culation is cumbersome although it is a straightforward 



task. The first term in Eq. (|A.ip includes the mean part 
go. It is easily obtained by substituting a — A. The 
fluctuation correction is obtained from the integral 



2 r°° 

Tg 



ds\ \ dse 

i J-i 



x 



Ti 



[l + f( Sl -l)] [l-f(l- S )] 
1-f [1 + f + 



;i + f( Sl -l)] [1-^(1 -a)] 



(A.8) 



In the same way, the second term in Eq. (1A.1|) is reduced 
to 



Tel dsi L ds J^sr' 

Ti 



To 



1 + #(«!-!) l + f( Sl -l) 



2 



T 



1-^(1- s) 1 + ^(1- s ) 



(l-s 2 ) 



(A.9) 



We note that these expressions are obtained after inte- 
grating the Grassmann variables and changing the vari- 
ables as s± = coshes and s = cos# F . A careful manip- 
ulation is required to carry out the remaining integrals. 
After lengthy calculations we can obtain Eq. (j34|) . 

It is a straightforward task to include the dephasing ef- 
fect described by the dephasing term |5T|) . In the present 
parametrization, it can be written as 



F r = T-( s l - S), 

At 



(A.10) 



and is incorporated in the integrals in Eqs. (|A.8p and 
(|A.9[) as e~ Fr . Although we do not show the analyti- 
cal result explicitly, it is not difficult to carry out the 
integrals. At the limit \E — Ei\ — > 00, we can find the 
result of Ref.fll. 
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